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selected and solved in a manner that builds up the use of 
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three sample problems to illustrate further the programming 
techniques presented in Section I« These problems 
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necessary computer program. The second is a problem frcm 
advanced algebra. The last is from the area cf statistics 
and represents a useful program for students (or teachers) 
for analyzing data. Related exercises fcr the reader are 
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duction 



C OMPUTERS have great promise as a tool in the teaching of mathe- 
matics. Many applications of the computer as a tool for learning 
through problem solving have been found, and many programs in com- 
puter education are in existence throughout the United States. It can be 
anticipated that the use of computers will become an integral part of the 
secondary school mathematics curriculum. Indeed, such is already the 
case in many school systems. 

This booklet aims to help the mathematics teacher introduce com- 
puters through an easy, problem-oriented language. In Section I a 
problem is selected and solved in a manner that builds up the use of the 
language. The statements of the language are few in number and can 
be clearly understood and readily learned. Then in Section II the lan- 
guage is applied to three additional mathematics problems. 

All the problems in the booklet can be solved on any computer that 
can accept the language called “BASIC.” It is expected that more and 
more computer facilities will make BASIC available to users. For facili- 
ties that do not accept this language, transfer to other usable languages 
is not difficult. After students are introduced to concepts of pro- 
gramming and can produce logical solutions to problems, the language 
becomes secondary. The knowledge transfer to other languages is very 
high. 




ducing BASIC 



I N INTRODUCING the programming language BASIC, a reasonable 
, approach is to consider an elementary but interesting mathematical 
setting. One interesting problem for the application of computers is to 
develop a computer algorithm for finding a good approximation to 
Vn for n > 0 — in computer typography, to Vn for n > = 0. 

Introductory Problem 

Normally, a satisfactory approach to a problem such as this is to 
consider first a special case for n, say 7. We will write a computer pro- 
gram to enable us to use the computer to find a good approximation to 

V?. 

A note describing a computer program is appropriate here. A 
program is simply an ordered sequence of statements or steps. For 
example, if you were to write out the rules of the traditional steps for 
long division, this could be thought of as a program even though it is 
not a formal computer program. 

In writing computer programs it is necessary to use a language or 
vocabulary that a particular computer “understands.” There are many 
computer languages. In this booklet we have selected one language, 
BASIC, over the others because of the small number of commands needed, 
its easy application in solving problems, and its practicality in our evolv- 
ing educational setting. At the present time there are many languages 
with much more extensive commercial use than BASIC. For example, 
FORTRAN and ALGOL are such languages. 

Now back to our initial problem — how to use the computer to find 
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V7 Our first programs will involve very few programming statements 
and hence will be Ly to write. Then, since these 
efficient in their use of the machine, we 

the elementary BASIC language statements to write . 

The improved programs will introduce you to additional statements 
BASIC and the sequence of program development should suggest a 
fediffique appmpril for teaching other mathematical topics with a 
commrter The ability to develop a general algorithm is, m fact, one of 
the broad goals of mathematics education and is considered by many 

f. n i-, e flag highest level of problem solving. . 

t0 before approaching the given problem a brief review note is ne , 
saiy A positive number x is the square root of n if x- = N. Using t 
notion we cm start with our first computer program. Initially we shall 
se the computer as a computational device, similar to a calculator. In 
the language BASIC the computer can he instructed to do compu a ion 
S Memcnt called a mnt statement A sample program will 
S k, A preliminary estimate by mental arithmetic , nd.cates that 
V7 is between 2 and 3 and is in fact greater than 2.5, since (2.5) — ^ 
Thus consider 2.6 as an approximation to V?. One way to check die- 
curacy of this approximation is to square 2.6 The foUowing BASIC 
nroeram tells the computer to do this computation and print the resul . 



1 Program A 

1 print 2 6 * 2.6 This command, or statement, instructs the com- 

puter to compute and then print the result of 
multiplying 2.6 by 2.6. (Note that the symbol ” 
is used for multiplication and stands for “times, 
so that 2.6 * 2.6 means “2.6 times 2.6/’) The 
symbols for addition and subtraction are the 
usual ones, namely, + ailc ^ A slash, / , is used 

for division. The notation for an exponent is a 
vertical arrow; for example, 2.6 t 2 means “2.6 
squared.” Thus Statement 1 could be written, 
alternatively, as “1 print 2.6 t 2.” 

0 END This statement tells the machine that this is all 

there is to this program. All programs in BASIC 
must end with an end statement. 

After a program has been read into the computer 
the operator types run when he wishes the com- 
puter to execute that program. 

One further note about our first program is ill order. You will have 
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noticed that each program statement, up to and including end, has a 
number in front of it. These numbers tell the computer the order in 
which the statements are to be executed— and this is always done from 
smallest to largest. The end statement must always have the largest pro- 
gram number. The numbers in our program did not have to be 1 and 2; 
t iey mig it bettei have been 10 and 20. Using these numbers would have 
enabled us to insert steps in between, if necessary, merely by calling the 
new steps 15, 17, etc. The steps can be typed in any order in the program, 
since the computer will always execute in numerical order; that is, the 
machine will automatically order the program statements by number, 
lins is a characteristic of many computer languages. 

Now let us consider the output or “printout” from the computer for 
our first program. This is given below. ( Each program written in this 
section will be repeated when the computer output is given for your 
information. In addition, the programming statements introduced up to 
that point will be listed. ) 

Program A 

1 print 2.6 * 2.6 

2 end 
RUN 



Output 

6.76 

BASIC programming statements introduced thus far: print . . . and end 

The approximation of 2.6 appears to be a reasonable guess, but is 
this as good an approximation as we can get? We can repeat the proce- 
dure for 2.7 or 2.8, if we wish. What about approximating to hundredths? 
Could we write the program to do some further calculations? One handy 
feature of the print . . . statement is the fact that many different expres- 
sions can be evaluated within the same statement, as long as they are sepa- 
rated by commas. Thus our next program could look like the following. 

Program B, 



5 print 2.60 t 2, 2.61 t 2, 2.62 t 2, 2.63 t 2, 2.64 t 2 2 65 t 

6 print 2.66 t 2, 2.67 t 2, 2.68 t 2, 2.69 t 2 2 70 t ’ 

10 END ’ 



Note that the two Piu.vr . . . statements are needed, since each line 
must be preceded by a command. 
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The program is repeated below, and the output is given. 

Program B, 

5 print 2.60 1' 2, 2.61 f 2, 2.62 t 2, 2.63 f 2, 

6 print 2.66 f 2, 2.67 f 2, 2.68 t 2, 2.69 f 2, 

10 END 

RUN 

Output 

6.76 6.8121 6.8644 6.9169 6.9696 7.0225 

7.0756 7.1289 7.1824 7.2361 7.29 

BASIC programming statements introduced thus far: print . . . (with 
commas ) and end 

Note that we have to refer back to our list of calculations in the 
given program to determine which x and x- go together. Also, if we were 
to repeat this type of calculation for a large number of approximations, 
it would get pretty tiresome to keep writing t 2. Since we are squaring 
every time, it would be handy to have some way to read a value from a 
list and then have a single command that instructs the computer to square 
the number. In BASIC we have a pair of commands that enable us to 
read data from a list. These are appropriately called the read . . . and 
data . . . statements. Read the following program and discussion carefully. 



2.64 f 2, 2.65 f 2 
2.70 t 2 



Program B ^ 

10 read x This instructs the computer to find a data list 

and set x equal to the first value on the list — 
in this case, 2.60 would be the first number in the 
data list. In BASIC any letter or any letter and 
a single digit is acceptable for the name of a 
variable (hence we could have used s, si, q, q2, 
or any other such symbol for a variable). 

20 data 2.60, 2.61, 2.62, 2.63, 2.64, 2.65, 2.66, 2.67, 2.68, 2.69, 2.70 

The data statement can fall any place in the pro- 
gram. This is the list of possible replacements 
for x. The computer goes to this list only when 
instructed to do so in a read . . . statement. Other 
wise, the computer pays no attention to the state- 
ment. For each successive read . . . command the 
computer selects the next piece of data from the 
list until out of data, and then the computer 
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automatically goes to end — this is a distinct char- 
acteristic of BASIC. (Note: The data must be in 
decimal or integer form; fractions and radicals 
are not accepted . ) 

30 print x t 2 This instructs the computer to print the value of 

x squared. 

40 go to 10 At this point it is necessary either to repeat the 

earlier read . . . and print . . . instructions or, if 
possible, go back in the program and use the 
earlier statements over again. We use a go to . . . 
statement to jump from one statement to another 
in a program. In our program the go to 10 sends 
the computer back to Statement 10, the read 
instruction. The computer repeats this process 
until out of data in 20 and then goes to end - (the 
highest numbered statement). 

50 END 

The above program is repeated below, and the output is given. 

Program B z 

i 1 ' 

10 READ X 

20 DATA 2.60, 2.61, 2.62, 2.63, 2.64, 2.65, 2.66, 2.67, 2.68, 2.69, 2.70 

30 PRINT X t 2 

40 GO TO 10 

50 END 

RUN 

Output 

6.76 

6.8121 

6.8644 

6.9169 

6.9696 

7.0225 

7.0756 

7.1289 

7.1824 

7.2361 

7.29 

OUT OF DATA IN 10 









' HUM ■ ■ iww ■ Ml I I H l . g l - up 



• i . ' »«> ppippyisp 



6 / Introduction to an Algorithmic Language 

BASIC programming statements introduced thus far: print • • • > end, 

READ . . . , DATA . . . , GO TO . . . 

There may be some question as to whether program B, is really 
better than B x . To illustrate the superiority of the type of program illus- 
trated by Bo, we need only consider the case of evaluating a polynomial 
such as 3x 3 + x 2 — x for successive values of x, say 2.60, 2.61, . . . , 2.70. 
The print . . . statement for this problem would be a long monster. Try it! 
This illustrates the necessity for using variables in our programs. 

At this point let us stop and see what we can use the computer to do. 
We can actually write a program to instruct the computer to do any 
particular computation for a given set of data (as long as we can use an 
algebraic form to write the computation to be performed ) . This is a 
useful technique; however, there may be times when we will not want 
to print the results of all the computations. For example, in program B 2 
the computations 2.66 2 , 2.67 2 , . . . , 2.70 2 are all unnecessary, since the 
-y/7 is bounded between 2.64 and 2.65. It would be useful to be able to 
have the computer do a calculation for a given value for x and then stop 
and let us select a second approximation on the basis of the first results. 
For example, if you wished to find Vl,433 to the nearest tenth, you 
might try 35.0, and then on the basis of the computation of 35.0 2 you 
would try a value larger or smaller. It would be a lot of useless work 
to do the calculations of 30.0 2 up to 40.0 2 (or some other guessed upper 
and lower bounds), stepping by tenths. There is a BASIC statement that 
allows the programmer (you) to decide what value you want to substi- 
tute for the variable. This is the input . . . command. The following 
program shows how the input . . . command can be used in finding a/7. 

Program B., 

10 INPUT X 
20 PRINT X t 2 
30 GO TO 10 
40 END 

After reading (and storing internally) the program, the computer 
stops and waits for the programmer or operator to put in a value for x. 
After it receives this value for x, it prints the result of squaring x. And 
then what do you think happens? Note the use of the go to . . . statement. 
Without the go to . . . command in this program the computer would 
stop and the program would be completed. Then to do a second ap- 
proximation, the programmer would have to start at the beginning of 
the program again. In this case, however, the program instructs the 
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computer to go back to Statement 10. The machine will return to this 
instruction and stop and wait for the operator to select another value 
for x. The following output shows a running of this program. (The values 
to be substituted were arbitrarily selected. ) 

Program B, , 

10 input x 
20 PRINT X t 2 
30 GO TO 10 
40 END 

RUN 

Output 

? 2.6 
6.76 
? 2.7 
7.29 
? 2.65 
7.0225 
? 2.64 
6.9696 
? 2.645 
6.996025 

? STOP 

BASIC programming statements introduced thus far: print . . . , end, 

READ . . . , DATA . . . , GO TO ... , INPUT . . . 

Notice that we needed to use an additional comment, stop, in the 
actual running of our program. This particular program might run 
“forever” unless we tell the computer to stop. We can type one of two 
possible words to tell the machine that we have done enough work on 
this problem. The word used in the sample running of the program was 
stop. After stop it is possible to start the computer on a new program. 
The reader should refer to a BASIC manual to familiarize himself with the 
actual operation of the teletype and the input capabilities and commands 
of the facility. The words run and stop do not appear in the computer 
program but, rather, are used to tell the computer to start or stop. These 
are a part of the actual running of the equipment, but they are not con- 
sidered a part of the program. 

There is one additional advantage in our last program, B«, over the 
previous ones. Do you see what it is? We also checked 2.645 2 . We were 



The computer typed the question marks. We 
typed the values of x. The computer printed the 
values of x 2 . 
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not able to do this in any of our previous programs as they were written. 
This would be possible if we revised them to include this particular value. 
The input . . . statement enables us to make decisions on the basis of the 
most recent information and adjust our approximations accordingly. We 
are not bound by preconceived ideas. 




We now have a very good computer program for approximating 
the square root of any given number n except that the program, while it 
enables us to get a good approximation, is quite inefficient. The investi- 
gator must sit at the console and actually communicate with the computer. 
This presents problems if the computer facility is not immediately ac- 
cessible. Also, with this last program, the time spent on-line with the 
computer is excessive. Actually, the computer could do hundreds of 
calculations just during the time you spend typing in numbers. It would 
be valuable if we could have the computer check how close a particular 
approximation is and adjust the next trial value accordingly. To do this 
we need to have a statement that is referred to as a “decision” command. 
In addition, we want the machine then to change the trial value for x by 
some predetermined increment. The following program makes use of 
these ideas and introduces two more elementary BASIC statements. 

Program C 

10 let x = 2.0 The let . . . command enables us to set some 

variable x at a predetermined value. In this case 
we arbitrarily selected 2.0. 
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20 ifx t 2> 7 then 50 The if . . . then . . . statement is the “decision” 

statement. In BASIC the computer action at this 
point is just what is implied. The computer com- 
pares the value of x 2 with 7. If, as in our first 



30 let x == x + .1 
40 go to 20 



50 print x — .1, (x — 



case when x == 2, the value of x squared is not 
greater than 7 ( that is, it is less than or equal to 
7), the computer proceeds with the next step in 
the program, 30. If, on the other hand, x squared 
is greater than 7, then the computer is told to go 
to a particular statement number, in this case 50. 
Thus the computer is asked to check a given 
condition. If the condition is false, the computer 
continues with the program; and if the condition 
is true, then the computer goes to the statement 
number given in the command. A further note 
is necessary: The condition to be checked can in- 
volve any one of the symbols that are shown in 
the box on page 11, together with those intro- 
duced earlier. It is up to the programmer to 
decide which condition is of interest. 

Statements 30 and 40 illustrate a very interesting 
technique. The computer in this situation is told 
to replace x by the value x + .1. Since x is 2.0 
prior to Statement 30, x is changed to 2.1, and 
then Statement 40 sends the computer back to 
Statement 20 to repeat the comparison of x 2 with 
7. Note that in Statement 30 we do not have a 
true application of the concept of equality. 
Rather, the equality symbol can be thought of as 
being a symbol for the notion “is replaced by.” 
This is a technique that is used in all of the 
computer languages. The technique is often re- 
ferred to as a ‘counter.” Note that if the in- 
crement on x were 1, and x started at 1, this 
would in fact make x equal to 2, 3, 4, and so on, 
for each “loop” of the program. (A section of a 
program that is repeated a number of times is 
referred to as a loop.) 

.1) t 2, x, x t 2 

When x squared first becomes greater than 7 the 
condition x 2 > 7 in Statement 20 is true, and thus 
the computer goes to Statement 50 given in the 
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then . , . part of the condition. In Statement 50 
we might have decided to print only x 2 , but ac- 
tually we do not know what the corresponding 
value of x is at this point ( all the computations 
and comparisons have been done in the machine). 
Therefore it is desirable to print x also. But what 
is true about x and x 2 ? They are the first values 
that are too large. Thus it would be useful, for 
this program, to print also the last value of x 
which still satisfied the condition jtljat x 2 < 7. 
Then we can decide which approximation to use 
for a specific purpose and whether to seek other 
approximations. Program procedures such as this 
are actually up to the ingenuity of the pro- 
grammer. Care should be taken to be sure that 
the computer output is in a readily applicable 
and usable form. In this case notice that we have 
bounded the true value of ■s/l to the nearest 
tenth. 

60 END 

The above program is given again, below, and the actual computer 
output is shown. The actual running time of the computer for this pro- 
gram is such that the final answer is typed almost immediately after the 
program has been put into the machine. The calculations are all done 
in less than half a second. This program is considerably more efficient 
than the previous program, B ;i . 

Program C 

10 let x = 2.0 
20 if x t 2 > 7 then 50 
30 LET X = X -f- .1 
40 go to 20 

50 PRINT X — .1, (x — .1) f 2, X t 2 

60 END 

RUN 

Output 

2.6 6.76 2.7 7.29 

BASIC programming statements introduced thus far: print . . . , end, 

READ . . . , DATA . . . , GO TO . . . , INPUT . . . , LET . . . , IF . . . THEN . . . 
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Meaning of Computer Symbols 

+ plus 
— minus 
* multiply by 
/ divide by 
t exponent 
< less than 
> greater than 

< = less than or equal to 

> = greater than or equal to 
= equal to ( or, is replaced by ) 

< > not equal to ( this symbol is used 

rather than ^ or =/) 

We might wish to stop now — or maybe you have an idea for a further 
refinement of the program for this problem. What if we wanted an 
approximation to the nearest .0001? Would you start at, say, 2.6000 and 
increment by .0001? How many calculations would be involved? (In this 
case there would be 458. ) What if we weren’t sure about 2.6 and started 
at 2.0 as in our previous program? — Then 6,458 calculations would be 
necessary. Thus the problem facing us is as follows: Can we reduce this 
number of computations and comparisons in some reasonable way and 
still have the computer make all the computations and decisions neces- 
sary? How about changing the increment? 

Study the following program carefully. This program illustrates one 
possible procedure for getting an extremely accurate approximation with 
a small number of calculations. 

Program D t 

10 LET X = 2.0 

20 LET K = .1 

30 IF X t 2 > 7 THEN 55 

40 LET X = X + K 

50 go to 30 

55 IF K < .001 THEN 90 

60 LET X = X — K 

70 letk = k/10 
80 go to 30 

90 PRINT X — K, (x — k) t 2, X, X t 2 
100 END 
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2.6457 6.99973 2.6458 7.00026 

Do you see how this program gets a closer and closer approxi- 
mation to y/f? Note that we have let k be our increment; as soon as 
(x 4" K )" becomes greater than 7, we drop back to the value of x just 
prior to the last incrementing and increase this value now by k/10. Thus 
we have at most ten calculations for each of the increments of .1, .01, .001, 
,0001. When k becomes .0001, what happens? We have now bounded 
the approximation to the nearest ten-thousandth; and, rather than go 
further, our decision statement (55) instructs the computer to go to 
Statement 90, the print statement. Read through the program carefully 
until you understand all parts, for this program represents as good a 
refinement of our original problem as we will attempt in this booklet. No 
new statements are introduced in Program D l} but it is somewhat longer 
and more complicated than any previously attempted. 

At this point in our discussion we have noted all but two of the ele- 
mentary BASIC statements. Addition of these two new statements, the 
for . . . and next , . . statements, allows us to replace many of the state- 
ments in the previous programs. The new statements set up the “counter” 
referred to in our discussion of loops. The for . . . statement enables us 
to set a variable equal to some value and increment it for each successive 
loop in the computer program. Consider the program statement given 
below: 

5 for x = 0 to 48 step 2 

This statement, in which step means to add an amount other than 1 (in 
this case, 2) should be read as though there were a comma separating 
the “48” and the “step.” It says to start with x = 0, proceed with the rest 
of the program, and increment by 2 each time the computer returns to 
this statement until x > 48. At some point in the program there needs to 
be a statement like the following: 

25 next x 

This tells the computer to return to the for . . . statement, increment ( that 
is, increase x by 2 as specified), and continue with the program again. 
Thus in Statement 5 the next x value will be 2, then 4, and so on. The 
machine will automatically check each statement and go to the statement 
that follows the next x statement as soon as x > 48. We can change the 
increment to whatever we wish. For example, Statement 5 could be as 
follows: 

5 FOR X “as 0 to 48 step .5 
or 
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5 for x = 0 to 48 step 7.25 

If we wish to increment by 1, then the step part of the statement can be 
left off, and the statement becomes 

5 for x — 0 to 48 

Remember that each for . . . statement in a program must have an ac- 
companying next. ... It is also possible to write the for . . . statement 
with variables — say a, b, and c: — and then let a, b, and c take on selected 
values. In this case the for . . . statement would look like the following: 

5 FOR X ss= A TO B STEP 0 

The two new for . . . and next . . . statements can be used to replace por- 
tions of Programs C and D t . An example follows, after we repeat Pro- 
gram C for the sake of ready comparison. 

Program C 

10 let x — 2.0 
20 if x t 2 > 7 then 50 
30 let x — x + .1 
40 go to 20 

50 print x — .1, (x — .1) t 2, x, x t 2 

60 END 

RUN 

Output 

2.6 6.76 2.7 7.29 

Program D,< 

( Program C with for . . . , next . . . ) 

10 FOR X ass 2.0 TO 3.0 STEP .1 
20 IF X t 2 > 7 THEN 40 
30 NEXT X 

40 PRINT x — .1, (x — .1) t 2, X, X t 2 
50 END 

RUN 

Output 

2.6 6.76 2.7 7.29 

Note that in Program Do we did not want to process the entire range 
of values set up in the for . . . statement. Instead, we gave the computer a 
rule for determining whether or not it was useful to continue processing 
values by means of the for . . . statement. We did this by inserting an 
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if . . . then . . . statement between the for . . . statement and its associated 
next .... In this way, we “jumped out of the for . . . loop before it 
would ordinarily end, thereby saving computer time. It is important to 
note that although we could “jump out of” the for . . . loop, we must start 
at the beginning in order to reenter the loop. 



BASIC programming statements introduced thus far: 

READ . . . , DATA . . . , GO TO . . . , INPUT . . . , LET . . . , 
FOR . . . NEXT ... 



PRINT . . . , END, 
IF . . . THEN . . . , 



In this new computer program notice how we use the for . . . and 
next . . . statements to replace three statements in the original program 
Program D lf reprinted below from page 11, offers even more of a chal- 
lenge— you might wish to try to write it, using for . . . and next . . . state- 
ments, before going on to the next section, so that you can compare your 
results with Program D :t , a sample program using for . . . and next. 

Program D, 

10 let x = 2.0 

20 let k = .1 

30 IF X t 2 > 7 THEN 55 

40 LET X = X -f K 

50 go to 30 

55 IF K < .001 THEN 90 

60 LET X = X — K 

70 LET K = x/10 

80 go to 30 

90 PRINT X - X, (x - k) t 2, X, X t 2 
100 END 

RUN 

Output 

2.6457 6.99973 2.6458 7.00026 

Program D, 

( Program D* with for . . . , n/jxt . . .) 

10 let x = 2.0 
20 FOR J = 1 TO 4 
30 LET x = 10 t ( — j) 

40 FOR X = X TO 3 STEP X 
50 IF X t 2 > 7 THEN 70 
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60 next x 
70 LET X = x — K 
80 NEXT J 

90 PRINT X, X t 2, X + K, (x + k) t 2 
100 END 

RUN 

Output 

2.6457 6.99973 2.6458 7.00026 

Note that m Programs D t and D ;t we actually are telling the com- 
puter to do the same computations; however, use of the for . . . and 
next . . . statements often shortens a program. Look at the program that 
uses the for . . . and next. ... Do you see what is happening? First x 
is set to a given value, 2, and then k is set at .1 (Statements 10, 20, and 
30 m Program D ; , actually do this for k). Then x is incremented by .1 
2?* 1S Sweater than 7. At this point the program goes to Statement 
70, which sets x equal to the value just preceding ths point at which the 
squared quantity exceeded 7. Then the increment is changed from .1 to 
.01 (this is actually 10'J, which is now 10' 2 ). Then the squaring and com- 
paring process is again done until x 2 is greater than 7 and we jump from 
the loop on x to die “outside” loop on x. These two loops are referred to 
as nested loops, which actually means a loop within a loop. 

With the addition of these two new statements it is possible to 

program many of the problems or algorithms in the secondary mathe- 
matics curriculum. 

A list of the elementary BASIC statements and the program in which 
each was introduced is given in the table below. Some of the statements 
are listed together, since they must both appear in the program. 

Sample Program 

A, Bi 
A 
B 2 
Bo 
B:t 

C 
C 

Do, D 3 

The letters A, B, C, and D are used to name the sample programs to 

indicate levels of refinement, the D ? s being the most sophisticated or 
best programs. 



Statements) 
print . . . 

END 

read . . . and data . . . 

GO TO . . . 

input . . . 

LET . . . 

IF . . . THEN . . . 

for . . . and next . . . 



1 
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Most computers have a stored program foi finding square root. ^ In 
BASIC this square root function, generally referred to as a subroutine, 
can be called in with the letters sqr followed by parentheses enclosing 
the number to be operated on. Thus we can eliminate all the previous 
programming and merely write the following piogiam. e ou pi 
also given, for your information. 



Program D, t 

10 PRINT SQR(7) 
20 END 

RUN 



Output 

2.64575 

Exercises 

(Sample programs for all exercises in Sections 1 and 11 are given in 
the concluding section of this booklet.) 

1. Write a program to extend C, Di, or Du so that the computei will 
read any value for n from a data list and then find the square root by 

successive approximations. , 

2. Check your results in 1 above by using a read . . . data . . . combina- 
tion and the square root function. 

3 Often we use an iterative process to find the square root of a number. 
One iterative process is to take your first approximation for the square 
root of n, call this x, and divide n by this approximation Call the 
quotient q. Then average q and x, take their average as the second 
approximation, and repeat the process. Write a computer program to 
do this to find V7 to four decimal places. 




e Problems 
xercises 



T HE following section contains three sample problems to illustrate 
further the programming techniques presented in Section I. These 
problems illustrate three different mathematical settings that lend them- 
selves to computer analysis. The first is from number theory and forces 
the student to recall precise definitions of certain key mathematical con- 
cepts in order to write the necessary computer program. The second 
is a problem from advanced algebra. The last is from the area of statistics 
and represents a useful program for students (or teachers) for analyzing 
data. The three problems are independent of one another and can be 
read in any order. At the end of each sample problem there are a few 
related exercises for you to try. Answers to all will be found in the back 
of this booklet. 

Problem 1 — Divisors of a Positive Integer 

Now that we know something about the BASIC language, we will 
use it to help us solve the following problem: 

Let m be a positive integer. Compute the set of positive integral 
divisors of m. 

Lets look at a specific value of m — say 36. The divisors of 36 are 1, 
2, 3, 4, 6, 9, 12, 18, and 36. We can obtain these divisors in a very short 
time by using paper and pencil. But how would we find the divisors 
of 3,942? One way is to divide 3,942 by 1, 2, 3, 4, 5, and so on, each time 
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examining the remainder. If the remainder is zero, then we have found 
a divisor. 

We will develop a BASIC program directing the computer to do 
most of the work for us. To do this, we must first select variables to rep- 
resent the various results of our divisions: that is, we wish to represent 
the integral part of a quotient and also the remainder. For any integral 
trial divisor n > 0 let us consider Q to be the integral part of m/n and 
r to be the remainder, where 

m , R 

Q+ V 

Therefore, for any positive integer M, consider N to be any consecutive 
integer from 1 to m — that is, 1, 2, 3, . . . , m. Since division can be checked 
by multiplication, the dividend can be expressed further in terms of n, 
q, and r — specifically, m = n • Q + R. Thus r = m — n • Q. To clarify this 
notation further, consider the following example: If m = 48 and N = 5, 
then q = integral part of 48/5 = 9; r = 48 — 45 = 3. 

Does m = n • q -J- r? 

Check: 48 = 5 • 9 + 3. 

In this case h 7^ 0; therefore 5 is not a divisor of 48. However, if N = 6, 
then q = integral part of 48/6 = 8, and r = 48 — 6 • 8 = 0. In this case 
r = 0; therefore 6 is a divisor of 48. 

One approach to obtaining the set of divisors of m is to write a 
BASIC program directing the computer to divide an input value of m by N 
for n = 1, 2, 3, ... , m. After each division we will have the computer 
print the value of n and the value of the remainder r. Hence we can look 
at the output and determine which values of n are divisors of m. 

In order to obtain the required integral part of the quotient on divi- 
sion of m by n, we will make use of a special function that is available to 
us. This is a stored program, as was the square root function, sqr( ), 
described in the previous section. This is the integer function, which is 
part of the BASIC language. This function can be elicited with the 
notation int( ) where the expression within the parentheses is to be 
operated on. If e is any algebraic expression, then int( e ) = the great- 
est integer < e. For example, 

int(5) = 5. int(5.8) =5. int( 3.6 -j- 2.7) = 6. 

int( — 6) = —6. int( — 6.3) = — 7. int(48/5)=9. 

int(0) =0. int( 57.392) = 57. int(48/6) =8. 

Hence int(m/n) is the integral part of the quotient on dividing m by n. 

Now, since we are able to represent q and r, we are ready to write 
a BASIC program to find the division of a positive integer m. For our 
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first program, let us consider m to be 10. The following program is one 
possible approach to this problem. 



Program E t 



10 INPUT M 

20 FOR N = 1 TO M 
30 LET Q = INt(m/n) 
40 LET R == M — N * Q 
50 PRINT N, R 
60 NEXT N 
65 GO TO 10 
99 END 

RUN 



We must first set our value for m. Then, using 
a for . . . next . . . loop, we can generate our set 
of trial divisors and compute the integer quotient 
q and the remainder r (Statements 30 and 40). 
The values of n and r are printed, and the com- 
puter is instructed to continue this process for 
all trial divisors. 



Looking at the output, can you see which are 
divisors of 10? The divisors of 10 are 1, 2, 5, and 
10, since in each of these cases r = 0. 



5 0 

6 4 

7 3 

8 2 

9 1 

10 0 
? STOP 



Output 

? 10 
1 0 
2 0 

3 1 

4 2 



Some of you may have been surprised by the computers ability to 
understand what we meant by the expression m — n * q in Line 40 of 
the program. After all, it is possible to interpret it as either “multiply q 
times the difference of m and n” or “subtract the product of n and q from 
m. How can we be sure that the computer will correctly interpret our 
intentions? In cases such as this the computer has a set of rules, which it 
will always follow. When parentheses are not used to indicate otherwise 
(and this, also, is part of the set of rules), it will perform all exponentia- 
tions first; next it will perform all multiplications and divisions, and finally 
it will do all additions and subtractions. If the total number of multipli- 
cations and divisions in a given expression is larger than one, then the 
computer will perform the leftmost multiplication or division (regardless 
of which it is ) first, skip over any additions or subtractions, and perform 
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the next multiplication or division from the left. When all multiplications 
and divisions are complete, this process is repeated for additions and 

subtractions. . , 

Hence in Line 40 the computer first looked for exponentiations, bind- 
ing none, it looked for a multiplication or division. It found n * Q and 
executed the multiplication. It found no further multiplications or divi- 
sions and hence looked for additions or subtractions. It found m — and t le 
product it had already computed, and it performed that subtraction. 
Hence it found R to be the difference of m and the product of N times Q, 

as we desired. nr ^ Ari . 

This same program can be rerun for the case M = 394-. The only 

change necessary is in Statement 10, which will use 3942 for m. A partial 

output for this program is given below. 



Output 
? 3942 



1 0 

2 0 

3 0 

4 2 

5 2 

6 0 

7 1 

8 6 

9 0 

10 2 

11 4 

12 6 

13 3 

14 8 

15 12 

16 6 

17 15 

18 0 

19 9 

20 2 



Now we might wish to develop a more elegant procedure to generate 
the set of divisors of a positive integer. Our first procedure does not 
begin to tap the power of the computer! As in our square-root example, 
it would be nice to have the computer make some of our decisions. Since 
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the quantity of output can become quite excessive, as in the case for 
m = 3942, it would be more efficient to include only the values of n that 
are divisors of m in the output. For each such value of n, it would also be 
convenient to include the quotient q in the print . . . statement. Can you 
write the refined program? One possible attempt is given in Program 
below (36 is a completely arbitrary choice for m). 

Program E 4 

10 INPUT M 
20 FOR N = 1 TO M 
30 LET Q = INT ( M/N ) 

40 LET R = M — N * Q 
50 IF R < > 0 THEN 70 
60 PRINT N, Q 
70 NEXT N 
80 GO TO 10 
99 END 

RUN 

Output 
? 36 

1 36 

2 18 

3 12 

4 9 

6 6 

9 4 

12 3 

18 2 

36 1 

? STOP 

Something rather interesting occurred in the above output. Every 
divisor of 36 appeared twice! In fact, all the divisors were printed in the 
first six lines of output ( 1, 2, 3, 4, and 6 in the first column and 36, 18, 
12, 9, and 6 in the second). Why did this happen? Can we take ad- 
vantage of this to write a more efficient program? We can. 

First, we note (it can be shown) that if we print the divisor n and 
the quotient q for each trial value of n, then we need use only the 
following set of trial divisors: 

{ 1, 2, 3, . . . , INT ( Vm)}. 



If R =t^= 0, then the computer is instructed to go 
to Statement 70, and hence the print . . . state- 
ment is skipped. 
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That is, we can stop when we have reached the value of n equal to the 
integral part of the square root of m. ( It is interesting to note that stu- 
dents in the classroom quickly come to this generalization or one very 
similar to it. ) For large values of m, this reduction in output represents 
a tremendous savings in computer time (and cost!). We have incor- 
porated this in the procedure below. (Note the use of the read . . . 
data . . . statements rather than the input . . . ; this allows us to consider 
more than one replacement for m. In this case we have m equal to 36, 42, 
and 851, respectively. ) 

Program E., 

10 READ M 
15 PRINT M 

20 FOR N = 1 TO INT(SQR(M) ) 

30 LET Q = INT(m/n) 

35 LET R = M — N * Q 
40 IF R < > 0 THEN 60 
50 PRINT N, Q 
60 NEXT N 
70 GO TO 10 
80 data 36, 42, 851 
99 end 

RUN 



Output 



? 36 


1 


36 


The divisors of 36 are 1, 2, 3, 4, 6, 9, 12, 18, and 


2 


18 


36. 


3 


12 




4 


9 




6 


6 




? 42 


1 


42 


The divisors of 42 are 1, 2, 3, 6, 7, 14, 21, and 42. 


2 


21 




3 


14 




6 


7 




? 851 


1 


851 


The divisors of 851 are 1, 23, 37, and 851. 


23 


37 




? STOP 
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There are still many ways to make the procedure more efficient. You 
will discover some of these by working the exercises at the end of this 
section. 

Let us now extend the problem we have worked on thus far to con- 
sider a concept frequently encountered in the upper elementary or junior 
high school grades, that of finding the greatest common divisor ( G.C.D. ) 
of two positive integers. 

Let a and b be positive integers. How can we utilize the computer 
to obtain the G.C.D. of a and b? One way would be to print out the sets 
of divisors and look at them to determine the G.C.D. Before continuing 
let us establish our notation as follows: 

Let Da be the set of positive integral divisors of a and Db be the set 
of positive integral divisors of b. Then let c = Da 0 Db. The greatest 
common divisor of a and b is the largest element of c. For example, let 
a = 36 and b = 42. Then 

d a = {1,2,3,4,6,9,12, 18, 36}, 
d b = {1,2,3,6,7,14,21,42}, 
c = {1,2, 3,6}/ 
and G.C.D. =• 6. 

Now, since we already have a program to generate the set of divisors 
of a positive integer, let us use it to get the G.C.D. of any positive integer 
a and b. Here’s how: 

1. Use the program to generate the set of divisors of a (Set d a ). 

2. Use the program to generate the set of divisors of b (Set Db). 

3. By inspection, write down the set c = Da fl Db. 

4. Select the largest element of c. This is the G.C.D. of a and b. 

Below is the program for a = 910 and b = 798. 

Program F t 
10 INPUT M 

20 FOR N = 1 TO INT(sQR(m) ) 

30 LET Q = INT(m/n) 

35 LET R = M — N * Q 
40 IF R < > 0 THEN 60 
50 PRINT N, Q 
60 NEXT N 
70 GO TO 10 
99 END 

RUN 



O 
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Output 


? 910 


r 

i 


1 


910 


2 


455 


5 


182 


7 


130 


10 


91 


13 


70 


14 


65 


26 


35 


? 798 


1 


798 


2 


399 


3 


266 


6 


133 


7 


114 


14 


57 


19 


42 


21 


38 



The value of b was given to the computer, and 
the list of divisors was computed and printed. 



? STOP 

After getting the above output, we need only select the common divisors 
from die sets. In our example the set of common divisors is 

1 1,2,7, 14} , 

‘I’® to we might ask ourselves: Is there a better way to do this prob- 
lem? Can we get the computer to search and make the ecisions so ° n 
the results of concern are printed? This problem is given for you as 

exercise. 

Exercises 

1 The nrogram above causes int(sqr(m) ) to be evaluated for each 
time through die loop. A square root computationisrelahv^hme- 
eonsuming compared to many other operations in BASIC. Kewiite the 

program so drat int(sqr(m) ) is dement 

and assigned to the variable, s. Then use s in me roa 

2 We have used these trial divisors: 

n = 1, 2, 3, ... , int(sqr(m) ). 

That is, our program terminates when n > int(sqr(m) ). n equiva- 
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lent condition is to terminate the procedure when n > q. Write a 
program in which this fact is used to terminate the program. 

3. Write a program to determine if an input value of m is a prime num- 
ber. We assume that the input value for m is an integer ^ 2. ( A num- 
ber, m, is a prime number if and only if it has exactly two distinct 
divisors, 1 and m. Hence 1 is not a prime number.) If m is prime, 
print 1 and m. If m is not prime, print 1, m, and n, where n is a 
divisor of m other than 1 or m. 

4. Write a program for computing the G.C.D. of two positive integers, 
a and b, that prints only those divisors of A and b. (Hint: Use read a, 
b and direct the computer to print n only if n is a divisor of both 
integers. ) 



Problem 2 — Maximizing an Area 

The following is an example of a well-known type of problem in- 
volving a maximum value of a simple function. 

A farmer has a barn 37.0 feet long with a very long single-strand 
barbed wire corral on one end of the barn. He also has 525.1 feet of 
barbed wire to build a new two-strand rectangular corral along the side 
of the barn and extending out to the side of the old corral, using the barn 
and one strand of wire from the old corral in one side of the new corral 
as shown in the figure. It is then necessary to add one strand of wire from 
a to e and to put two strands along each of the line segments, ab, bc, 
and dc. 
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What should the dimensions of the new corral be in order for the 
corral to have the greatest possible area? 

If ab is x feet long and bc is y feet long, then ae is y — 37.0 feet long; 
and, adding the length of one wire for ae and two wires for each other 
side, we get the equation for the length of barbed wire used in the new 
corral: 
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2x + 2y + 2x + y - 37.0 = 525.1, 



or 



y = (562.1 — 4x)/3. 

The area of the new rectangular corral is a = xy, or, substituting for y 
from above; 



A = XY, 

562.1 - 4x 
A = x 2 , 

562. lx - 4x 2 
A =— 3 

(This is a quadratic function whose graph is a parabola opening down.) 
The reader should study this carefully to be sure he understands the 
derivation of this function. 

The questions now become: 

For what value of x (that is, length of ab) is the area of the corral 
a maximum? From the graph, what is the x-coordinate of the vertex of 
the parabola? 

Let us consider a very simple program to solve this problem. (You 
may wish to try a first attempt on your own. ) A first approach might be 
to input a value for x and then have the computer print x and a. The 
following program illustrates this approach. Values of 50, 70, and 90 
were arbitrarily selected for the sample output. 



The student makes an initial guess for x. 
x — 4 * x t 2)/3 

The computer prints x, the field length, and the 
calculated area. 

The computer is instructed to return to Step 10 
above for a new guess by the student. 

The end of the program. However, this program 
never gets to the end, since Step 30 always 
branches to Step 10. 



70 6582.33 



Program G t 

10 INPUT X 
20 print x, (562.1 1 

30 GO TO 10 
50 END 

Output 
? 50 

50 6035 

? 70 
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? 90 

90 6063 

? STOP 

This, of course, is a slow process, since there is a necessary wait while 
the student chooses and types in data. It may also, unless the student is 
a very astute guesser, take a long time to produce a nearly correct answer. 

As in the earlier examples, the question is how we can improve this 
program. Let us look more closely at the problem. If x = 0, then the 
area is 0. If x = 1, then the area is (562.1 • 1 — 4 • l 2 )/3 — that is, 
558.1/3. It appears that as we choose larger values for x, the area gets 
larger. But what about x very large, say 1000? Then a is negative. Thus 
the area reaches a maximum value for some x; then becomes smaller until 
x = 562.1/4 and the area is again 0. To verify this pattern, the reader 
should refer to the output from Program Gi. 

For our second approach to the problem, we might have the com- 
puter write out x and the area for every integral value of x in the set 

{ x: 1 <x< 141}, 

since 141 > 562.1/4. 

Program G 2 

10 for x = 1 to 141 This statement produces a loop that will be re- 
peated 141 times, each time increasing x by 1. 
Thus x takes on the values 1, 2, 3, ... , 141. 

20 print x, (562.1 * x — 4 * x t 2)/3 

The computer calculates and prints x and the 
area when the length is x, for each x from x = 1 
to x = 141. Each for . . . must have a next. . . . 
30 next x Each time this step is encountered the computer 

returns to Step 10, where x is increased by 1, 
until the process has been completed with 
x = 141. On the 141st pass the program goes on 
to the statement following the next . . . com- 
mand: in this case, Statement 40. 

40 END 

Part of the output 
1 186.033 
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5 903.5 

6 1076.2 

7 1246.23 

8 1413.6 

9 1578.3 

10 1740.33 

11 1899.7 

12 2056.4 

13 2210.43 

14 2361.8 

15 2510.5 



65 6545.5 

66 6558.2 

67 6568.23 

68 6575.6 

69 6580.3 

70 6582.33 

71 6581.7 

72 6578.4 



STOP 



It is still not clear what value of x provides a maximum area, but the 
solution must, be between 69 and 71, since in this interval the value for 
the area reaches a value greater than that for either 69 or 71 and begins 
a continuing decrease. We can now write a program that would search 
between 69 and 71, with smaller intervals, for a better approximation of 
the correct x value. This would entail only a change in the for . . . state- 
ment. The object of this example is to design a better program to search 
for x. 

In the next paragraphs in this section we will set up a program that 
will determine for itself the interval for the finer search and then search 
it. Notice that this is quite similar to our task of narrowing down our 
estimate for V7. We can make this a continuing process so that each 
search of an interval yields a new, smaller, interval for a new search with 
even smaller subintervals. 
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Before considering our program, a word of caution about the prob- 
lem is in order. It might have been assumed in the preceding discus- 
sion that, because when x = 71 the area is less than the area when 
x = 70, the maximum area would be found somewhere in the interval 
70 < x < 71; this is not necessarily true, although it may be. 

We can assume that the area function looks something like the figure 
below, with s, t, and u consecutive integers. In this example, even 
though the maximum area is reached between the integers s and t, the 
first noted decrease is between the integers t and u. In order to be cer- 
tain that the chosen subinterval will include the maximum value for the 
area, it will be necessary that the subinterval be s < x < u. That is, 
the subinterval must be two units long, these units being the last interval 
where the area increases and the first interval where the area deci eases. 




With this information in mind, our new program should do the 
following: 

1. Search from x 1 toward x — 141 looking for the first decrease 
in area. 

2. Type out the end coordinates (x, area) of the interval on the 
curve that includes the maximum area, with the x interval two spaces 
long. 

3. With units one tenth as long as the units of the previous search, 
the program will search the new interval (end points given above) for 
a decrease in area and again type out coordinates. 

4. Return to Step 3 three times, each time finding and typing the 
end points of shorter and shorter subintervals containing the maximum 
value of the function. 

5. Stop. 

Referring to the graph given earlier, then: 

1. The interval of search will always be s < x < u. The first value 
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of s will be 1 and of u, 141; but these will change as new subintervals 
are chosen. 

2. r will always be the length of the subinterval into which we are 
breaking su for the search. On the first search r = 1; for the second, 
r = .1; on the third, .01, etc. 

The final program for this problem is given below, followed by a 
detailed discussion of each of the statements. 



Program G :t 

10 LET S =1 

20 FOR I = 1 TO 3 

30 LET R = 10/10 t I 

40 LET a = (562.1 * s - 4 * s t 2)/3 

50 let b = (562.1 * (s + r) — 4 * (s + r) t 2)/3 

60 IF B > A THEN 75 

70 GO TO 110 

75 LET C = A 

80 LET A = B 

90 LET S = S + R 

100 GO TO 50 

110 PRINT S — R, C, S -f- R, B 
120 LET S = S — R 
130 NEXT I 
140 END 

Discussion of Program G, j 

10 let s = 1 We will set s.= 1, since for the first search this 

interval su is 1 to 141. 

20 for i = 1 to 3 This statement forms a loop that ends at Step 

130. This loop will be traversed three times, the 
first time through i = 1, the second through 
i = 2, the third through i = 3. ( This marks the 
beginning of the search of an interval.) 

30 let h= 10/10 t i When i = 1 (first time through the loop) this is 

r = 10/10 1 or 1; when i = 2, r = 10/10 2 = .1; 
when i = 3, r = 10/10 3 = .01. 

40 let a = (562.1 ~ s — 4 * s t 2)/3 

The computer calculates the value for a at the 
left end of the first subinterval. Each time a new 
subinterval is defined for a search, the program 
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will return to this step, where the first (right- 
end) value of the area will be computed. At all 
other times when a left-end value for the area is 
needed the old right-end value will be used and 
only a new right-end value computed, in order 
to save computer time. 

50 let b = (562.1 * (s + r) - 4 * (s + r) t 2)/3 

The computer calculates the area at the right 
end of the first interval, and the program returns 
here (every time the area does not decrease) to 
compute a new right end. 

b is compared to a, and the computer then ex- 
ecutes Step 75 if b > a; otherwise Step 70 is 
executed. See Section I for other possibilities 
for the if . . . then . . . statements using >, =, 
etc. 

When b < a we have discovered a subinterval 
two units long containing the maximum value, 
and the computer is instructed to go to the print 
. . . statement. 

This is done to save a as c for printing in case 
it is needed as the left end of a two-unit interval, 
should the next interval show a decreased area. 
Statement 80 moves the interval under consid- 
eration over one unit to the right by moving the 
value at the right end to the left end, and s is 
increased for use in the calculation of a new 
right end in Step 50. 

The computer now returns to calculate a new 
right end. Notice that this does not affect the 
loop established in Step 70, since the next . . . 
statement is not encountered. 

110 PRINT S — R, C, S -f R, B 

The computer reaches the print . . . when b < a. 
The print . . . includes the end coordinates of an 
interval two units long containing the maximum 
area. Since the value of s is now the value of x 
at the left end of the one-unit-long subinterval 
where the first decrease was noted, the value of 
x at the left end of the two-unit-long interval will 
be one subinterval back, at s — r. 



60 IF B > A THEN 75 

70 GO TO 110 

75 LET C = A 

80 LET A = B 
90 LET S = S -f- R 

100 GO TO 50 
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120 let s = s — R 
130 next i 



140 END 




Output 


69 


6580.3 


70.2 


6582.42 


70.25 


6582.43 


71 


6581.7 


70.4 


6582.4 


70.27 


6582.43 



This statement resets the starting point of the 
new, shorter, interval for the next search. 

The computer returns to the beginning of the 
interval search. Since this is the end of a for . . . 
loop, the value of i is incremented on the return 
to Step 20, and therefore in Step 30 the length of 
the subinterval is divided by 10 for the next 
search. After the third pass the i is exhausted; 
the computer will not return to Step 20 but will 
go on to Statement 140, the statement following 
the next . . . statement. 

The end of the program. 



Exercises 

1. Rewrite the above program to find the maximum value of the function 



f(x) = — 4x 4 - 3x 3 + 2x - 1. 

2. Using a for . . . loop, write a program to evaluate the polynomial at 
the following values for x : x = 1, 2, 4, 10, 30. 



x 2 



f(x) = i + * + at +7u + zr + •• • + 



x l 



^20 

onf 



Problem 3— Mean and Standard Deviation 

Given a set of test scores, write a program to compute the mean and 
standard deviation of the data. 

Review: The defining formulas for the means, designated as x, and 
the standard deviation, s, are given below: 

N 

2 >, 

i = 1 Xi -f- Xo -f- Xr -f- , . . + x n 
X = 



(i) 



N 



N 
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N / N 

(X| — x) 2 | / (X! x) 2 

i= 1 1/ i = l 



(Note: When sample data are used to estimate s 2 for a larger popu- 
lation, the denominator used is n — 1 rather than n. ) 

Generally the following computational formula is used for s 2 , rather 
than the definitional form given in (ii) above. This computational for- 
mula can be readily derived from (ii) by expanding and simplifying. 




and s is the square root of this quantity. 

With formulas (i) and (iii) one can compute the mean and standard 
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deviation by summing the values of X[ and Xf and then substituting 
these summations in the appropriate formulas. 

To approach the problem stated initially, we now need to write a 
program to take any set of input data and perform the indicated calcula- 
tions. You may wish to make an attempt at solving this problem before 
reading the following discussion. Consider the scores for this problem 
to be the following: 71, 65, 70, 85, 90, 95, 61, 73, 81, 88, 75, 75, 91, 66, 69. 
(Your own program will probably deviate slightly from the sample that 
will follow. ) 

To input the data given in our problem, an appropriate technique is 
to use the read . . . and data . . . statements. Now you must decide what 
you wish the computer to do with each data value. In this problem it 
seems our first step is to sum the x’s and the x-’s; we shall call these sums 
m and s, respectively. To have the computer do the summing, it is neces- 
sary to start with the values of m and s each equal to zeio and then in- 
crease them by the appropriate X| and s j - values. At this time it is also 
appropriate to input n; this can be done with a read . . . , an input . . . , 
or a let . . . statement. (Can you write the program now?) 

Your program should look something like the following one. 

Program H 

10 le t n = 0 A counter n is established to enable the computer 

to keep track of the number of scores in the data 
statement. 

Statements 20 and 30 set the variables m and s, 
arbitrary choices for variables, equal to zero. 
Now we wish to take each piece of data (in the 
data . . . statement, 50) and do the calculations 
given in Statements 70 and 80. 

50 DATA 71, 65, 70, 85, 90, 95, 61, 73, 81, 88, 75, 75, 91, 66, 69, 9999 

(Note that 9999 is not one of our original test 
scores. — The use of this value or one like it, 
which is often referred to as a “flag,” will be 
pointed out in the discussion pf the next state- 
ment. ) 

60 if x > = 9999 then 100 

Statement 60 asks if the value of x has reached 
the large “extra” value in the data . . . statement. 
In the event x ^ 9999 then the computer will 
continue with the next statement in the program. 
If, however, x ^ 9999 then the computer will 



20 let m == 0 
30 let s = 0 
40 READ X 
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jump to Statement 100, which is the first of our 
two print . . . statements, The technique of in- 
cluding the 9999 or some other “flag” is essential 
to the program. If such a technique is not used, 
the computer will automatically go to end when 
out of data in 60; that is, when asked to read x 
and there is no value left for x, the computer 

will go to END. 

65 let N = n -|- 1 This increments the counter n by 1. 

70 let m = x + m This statement increases m by x for each succes- 

sive value of x; that is, the final m will ^e the 
sum of the x’s. 

80 let s = x t 2 -|- s Do you see what happens to s? It increases by 

x- for each successive value of x. Thus the final 
s will be equal to the sum of the x’s. 

90 go to 40 Statement 90 sends us back to 40, where the next 

value of x is read, and the interim statements are 
executed again until x = 9999 is encountered. 

100 PRINT M, N, “mean = ,” m/n 

Statement 100 tells the computer to print m, n, 
the phrase “mean = ,” and m/n. (This is the 
application of formula (i).) Note that we can 
include letters or words in our output merely by 
enclosing them in quotation marks in a print 
. . . statement. You may not wish to print out all 
this information; this is up to the discretion of 
the programmer. 

110 PRINT S, N, “s.D. = ,” SQR( (n * S — M t 2)/(n t 2) ) 

Statement 110 instructs the computer to print 
the values of s and n, “s.d. = ,” and the square 
root of the number within the parentheses. The 
computations come from formula (iii). (Note 
that this last calculation could be simplified 
somewhat by taking just the square root of the 
numerator and dividing this by n.) The BASIC 
subroutine for square root was introduced in Sec- 
tion I of this booklet. Recall that to use this 
function it is only necessary to use the letters sqr 
and insert the quantity to be operated upon 
within parentheses. 

The program is repeated below, and the output is given. 
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Program H 

10 let n = 0 
20 let m = 0 
30 LET S = 0 
40 READ X 

50 DATA 71, 65, 70, 85, 90, 95, 61, 73, 81, 88, 75, 75, 91, 66, 69, 9999 

60 if x > = 9999 then 100 

65 LET N = N -j- 1 

70 LET M = X -{- M 

80 LET S == X t 2 -j- S 

90 go to 40 

100 PRINT M, N, “MEAN = ,” m/n 

110 PRINT S, N, “s.D. — ,” SQR( (n * S — M t 2)/(n t 2) ) 

120 END 

Output 

1155 15 mean = 77 

90519 15 s.d. = 10.2762 

The above program can be easily altered for a different set of data 
by rewriting Statement 60. 

Exercises 

1. Write a program to be used with the program already written that 
will translate each raw score Xi into an appropriate z-score. The 
z-score translation changes the given distribution into a new distribu- 
tion with mean 50 and standard deviation 10. The translational for- 
mula is 



Use s.d. = 10.3. (Hint: In this problem you will need to print . . . 
after each translation.) 

2. Write a program to print the sum of sets of n consecutive odd inte- 
gers — that is, 1, 1 -f 3 = 4, 1 -J- 3 + 5 = 9, etc. — for n = 1 to 15. 





T HIS booklet is an introduction to the use of a problem-oriented com- 
puter language in the teaching of mathematics. For the sake of 
simplicity, no attempt has been made to utilize all the available features 
of the BASIC language. These features are available to the reader who 
wishes to refer to the manuals listed in the bibliography. 

Besides BASIC there are currently in existence several other problem- 
oriented computer languages. Examples of these are JOSS, TELCOMP, 
ITTRAN, CUPL, ESI, CAL, INSTRUCTRAN, FORTRAN, and ALGOL. 
PL1 is a new one— available to a limited extent but still under develop- 
ment. The COBOL language, which is in common use, is primarily for 
the business type of problems. 

There are also many special program languages. LISP is especially 
designed for handling lists of elements, and it is one of many such. COGO 
and STRESS are languages that have been especially designed for civil 
engineering problems. It is highly probable that we will see in the future 
many more especially designed computer languages to solve particular 
classes of problems. 

“Flow charting” can also be called a language. This language is, in 
essence, a systematic and reasonably complete method of outlining the 
method and steps in solving a problem. It employs a sequence of boxes, 
circles, and other such closed figures connected and interconnected in 
such a way that the problem solution is laid out by their use. This is 
accomplished by having each box contain a verbal or mathematical ex- 
pression of one or more steps of the problem. This scheme has the ad- 
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vantage that it does not in any way depend on the particular machine 
to be used or on any particular machine language. 

Many people think of flow charting as a necessary prerequisite to 
computer programming. It is certainly an aid in programming complex 
mathematical problems. The degree of complexity of the problem and 
the skill of the programmer help determine the extent to which flow 

charts are used in problem solution. 

These languages are at the very heart of what is now called com- 
puter science, which is crystallizing out as a separate discipline and 
being included in the college curriculum. Thus teachers now have an 
opportunity to receive more assistance from colleges where this disci- 
pline is being incorporated. It is expected that these curricula will soon 
include explicit courses for teachers. 




T HE following bibliography is not intended to be complete. Rather, 
it is given to suggest sources of information relevant to curriculum 

utilization of computer systems. 

Albrecht, Robert L., and Mara, Walter. Computer Math 
for High School Reading, Mass.: Addison-Wesley Publish- 
ing Co., scheduled for publication August 1968. 

Control Data Corporation. The Teacher-Student Approach 
to Computer Programming Concepts , Vols. I cind II. The 
Corporation, 1963. 

W H. Freeman and Company. Information (a “Scientific 
American” book). San Francisco: The Company, 1966. 

General Electric Company. BASIC Language, Reference 
Manual The Company, 1967. 

. Introduction to Programming in BASIC: An Elemen- 
tary Instruction Guide. The Company, 1967. 

Hoffman, Walter, et al “Computers for School Mathe- 
matics,” The Mathematics Teacher, LVIII (May 1965), 
393-401. 

Johnson, Donovan A., and Rising, Gerald R. Guidelines to 
Teaching Mathematics. Belmont, Calif.: Wadsworth Pub- 
lishing Co., 1967. 
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National Council of Teachers of Mathematics. Computer 
Facilities for Mathematics Instruction. Washington, D.C.: 
The Council, 1967. 

. Computer Oriented Mathematics: An Introduction for 

Teachers. Washington, D.C.: The Council, 1963. 

School Mathematics Study Group. Algorithms, Computa- 
tion and Mathematics. Stanford, Calif.: Stanford University, 
1965. 



wers to Exercises 



Section I 

Introductory Problem 



Exercise 1 

To extend C, the program should look like the following. Extensions for 
Dj and Do make use of similar changes. 



Program 

5 READ N 
10 LET X = 2.0 
20 IF X t 2 > N THEN 50 
30 LET x = X + .1 
40 go to 20 

50 PRINT X — .1, (x — .1) t 2, 

x, x t 2 

60 GO TO 5 
70 DATA 7, 14, 21 
80 END 

50 PRINT x — .1, (x — .1) t 2, X, X t 



RUN 








Output 


2.6 ' 


6.76 


2.7 


7.29 


3.7 


13.69 


3.8 


14.44 


4.5 


20.25 


; 4.6 


21.16 



OUT OF DATA IN 5 



Note that this program fails for n < 4. 
You may wish to change the let . . . 
statement to overcome this. Also, the 
increment of .1 is quite small, if the 
first n is something like 9980. You 
may wish to improve your program to 
handle this. 



2 
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Exercise 2 

Pro grain 

10 READ X 
20 PRINT SQR(x) 
30 GO TO 10 
40 DATA 7, 14, 21 
50 END 

RUN 

Output 

2.64575 

3.74166 

4.58258 

OUT OF DATA IN 10 



Exercise 3 

Program 

10 READ N 

20 LET X = 2.0 

30 LET Q = n/x 

40 IF X > Q THEN 70 

50 IF Q — X < .0001 THEN 100 

60 go to 80 

70 IF X — Q < .0001 THEN 100 
80 LET X = (q + x) /2 
90 go to 30 

100 PRINT X, Q, X * Q — N 
105 GO TO 10 
110 DATA 7, 14, 21 
120 END 

RUN 

Output 

2.64575 2.64575 

3.74166 3.74166 

4.58258 4.58257 



Note the need for both Statement 50 
and Statement 70. The difference be- 
tween Q and x may be great but nega- 
tive, and the approximation would be 
such that this would be the printed 
result even though incorrect. 



-7.45058 E— 9 
-1.49012 E— 8 
-2.98023 E— 8 



O 

ERIC 
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Section II 

Problem 1 — Divisors of a Positive Integer 

Exercise 1 



Program 

10 READ M 
15 PRINT M 

18 LET S =INT(SQR(M)) 
20 FOR N = 1 TO S 
30 LET Q = INT(n/n) 

35 LET R — M — N * Q 
40 IF R < > 0 THEN 60 
50 PRINT N, Q 
60 NEXT N 
70 GO TO 10 
80 DATA 36, 42, 851 
99 END 

RUN 

Output 

36 

1 36 

2 18 

3 12 

4 9 

6 6 

42 

1 42 

2 21 

3 14 

6 7 

851 

1 851 

23 37 

OUT OF DATA IN 10 

Exercise 2 

Program 

10 READ M 
15 PRINT M 
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20 let n = 1 
25 LET Q = INt(m/n) 
30 IF N > Q THEN 10 
35 LET R — M — N ! ' : Q 
40 IF R < > 0 THEN 50 
45 PRINT N, Q 
50 LET N = N + 1 
55 go to 25 
80 DATA 36, 42, 851 
99 END 

RUN 

Output 



36 
18 
12 

9 
6 

42 
21 
14 
7 

851 

37 

OUT OF DATA IN 10 

I 

I 

j Exercise 3 

Program 
10 READ M 

15 IF M <>int(m) then .10 

j 20 IF M < 2 THEN 10 

25 LET s = int.(sqr(m) ) 

30 FOR N = 2 TO S 
35 LET Q = INT(m/n) 

40 LET R = M — N * Q 
45 IF R = 0 THEN 65 
50 NEXT N 
55 PRINT 1, M 
60 GO TO 10 

i 65 PRINT 1, M, N 

! 70 GO TO 10 

80 DATA 3.7, 0, 1, 2, 3, 4, 15, 17 
99 END 



i 

i 



o 

ERIC 



36 

1 

2 

3 

4 
6 

42 

1 

2 

3 

6 

851 

1 

23 



RUN 
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Output 

1 2 

1 3 

1 4 2 

1 15 3 

1 17 

OUT OF DATA IN 10 

Exercise 4 
(3 alternatives) 

Exercise 4A 

Program 

10 READ A, B 
15 FOR N = 1 TO A 
20 let q1 = int(a/n) 

25 LET Rl = A — N * Ql 
30 LET q2 = int(b/n) 

35 LET r2 = B — N * q 2 
40 IF Rl < > 0 THEN 55 
45 IF r2 < > 0 THEN 55 
50 LET G = N 
55 NEXT N 
60 PRINT A, B, G 

65 GO TO 10 

70 DATA 45, 36, 34, 26, 50, 25, 97, 6, 
99 END 

RUN 



Output 



45 


36 


9 


34 


26 


2 


50 


25 


25 


97 


6 


1 


25 


50 


25 



OUT OF DATA IN 10 



Exercise 4B 

Program 

10 READ A, B 

15 FOR N = A TO 1 STEP — 1 
20 LET Ql = int(a/n) 

25 LET Rl = A — N ,: ' v q1 
30 LET q2 = int(b/n) 

35 let r2 — b — N q 2 
40 IF Rl < > 0 THEN 50 
45 IF r2 = 0 THEN 60 



Note that no output occurred for 
m == 3.7, 0, or 1, since these values 
violated the conditions for input. 



If both remainders are 0, then n is a 
divisor of both a and b. Since n in- 
creases by 1 each time, g will be the 
largest such divisor at the termination 
of the for . . . loop. 

;, 50 



Since n starts at a and then decreases 
by 1 each time through the loop, the 
first time that both Rl and r2 are 0 oq 
division by n the result will be the 
desired g.c.d. 
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50 NEXT N 
60 PRINT A, B, N 

70 GO TO 10 

80 DATA 45, 36, 34, 26, 50, 25, 97, 6, 25, 50 
99 END 

RUN 



Output 



45 


36 


9 


34 


26 


2 


50 


25 


25 


97 


6 


1 


25 


50 


25 



OUT OF DATA IN 10 



Exercise iC 

Program 

10 READ A, B 
15 LET M = A 
20 LET N = B 
25 LET Q = INT(m/n) 

30 LET R = M — N * Q 
35 IF R = 0 THEN 55 
40 LET M = N 
45 LET N = R 
50 GO TO 25 
55 PRINT A, B, N 
60 GO TO 10 

70 data 45, 36, 34, 26, 50, 25, 97, 6, 25, 50 
99 END 

RUN 



Output 



45 


36 


9 


34 


26 


2 


50 


25 


25 


97 


6 


1 


25 


50 


25 


OUT OF DATA IN 10 



Problem 2— Maximizing an Area 

Exercise 1 

Program 

10 let s = 0 
20 for i = 1 to 3 
30 LET R= 10/10 f I 
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40 let a = ( (— 4 s — 3) ,! s - s + 2) ■' s 1 



50 let b = ( (-4 * p — 3) * p * p + 2) * p - 1 

60 IF B > A THEN 75 

70 GO TO 100 

75 LET C = A 

80 LET A = B 

90 LET S = S + R 

95 go to 45 

100 PRINT S — R, C, S + R, B 
120 LET s = S — R 
130 NEXT I 
140 END 

RUN 

Output 



.36 - .487153 .38 - .488021 

Exercise 2 

Program 

10 READ X 

20 LET S = X + 1 

30 LET F = 1 

40 FOR Z = 2 TO 20 

50 LET F = F * Z 

60 LET s = s + X f z/f 

70 NEXT Z 

80 PRINT X, S 

90 GO TO 10 

100 data 1, 2, 4, 10, 30 

110 END 

RUN 



45 LET P = s + R 



The polynomial can be written in this 
factored synthetic substitution form 
(see Statement 40), and this reduces 
the total number of operations to be 
performed; an additional factor to 
consider is that the computer is rela- 
tively slow in working exponentials. 
We use Statement 45 to save calculat- 
ing time and computing time. We 
calculate p once and then use p in- 
stead of s + R in Statement 50. 



-1 0 
.3 - .5134 



1 -6 

.5 - .625 



Output 
1 2.71828 
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2 7.38906 

4 54.5981 

10 21991.5 

30 3.77068 e 11 

OUT OF DATA IN 10 



Problem 3— Mean and Standard Deviation 

Exercise 1 

Program 
10 read x 

20 print 10 * (x - 77)/10.3 + 50 
30 GO TO 10 

40 data 71, 65, 70, 85, 90, 95, 61, 73, 81, 88, 75, 75, 91, 66, 69 
50 END 

RUN 

Output 

44.1748 

38.3495 

43.2039 

57.767 

62.6214 

67.4757 

34.466 

46.1165 

53.8835 

60.6796 

48.0583 

48.0583 

63.5922 

39.3204 

41.2621 

OUT OF DATA IN 10 



Exercise 2 
(2 alternatives) 

Exercise 2 A 



Program 



10 let s = 0 

20 READ X _ _ „„ 

30 data 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 2o, 27, 29 



40 let s = s + x 
50 print s 
60 go to 20 



I 



I 
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70 END 

RUN 

Output 

1 

4 

9 

16 

25 

36 

49 

64 

81 

100 

121 

144 

169 

196 

225 

OUT OF DATA IN 20 

Exercise 2B 

Program 
10 let s = 0 

20 for x = 1 to 29 step 2 

30 LET S = S + X 

40 PRINT S 

50 next x 

60 END 

RUN 

Output 

1 

4 

9 

16 

25 

36 

49 

64 

81 

100 

121 

144 

169 

196 

225 



O 
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